FERMILAB-PUB-1 1-200-T 
CERN-PH-TH/201 1-109 



A Forward Branching Phase-Space Generator 



Walter T. Giele*, Gerben C. Stavenga^ 

Theoretical Physics Department 
Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 



Jan Winter* 

PH-TH Department, CERN, CH-1211 Geneva 23, Switzerland 



June 27, 2011 



Abstract 

We develop a forward branching phase-space generator for use in next-to-leading order parton 
level event generators. By performing 2^-3 branchings from a fixed jet phase-space point, all 
bremsstrahlung events contributing to the given jet configuration are generated. The resulting phase- 
space integration is three-dimensional irrespective of the considered jet multiplicity. In this first study, 
we use the forward branching phase-space generator to calculate in the leading-color approximation 
next-to-leading order corrections to fully differential gluonic jet configurations. 
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1 Introduction 

The evaluation of the virtual corrections to processes including a large number of jets has become straight- 
forward ifTmOl due to the advent of algorithmic implementations Il20] - E7l of generalized unitarity II281I29I 
using parametric integration techniques Il30l . The only required diagrammatic evaluations are partonic 
tree-level amplitudes. These diagrams can be efficiently evaluated through algorithmic recursion rela- 
tions ||3TH42|. As a consequence, all building blocks exist to construct a Next-to-Leading Order (NLO) 
event generator for the evaluation of high multiplicity jet observables. Such a generator would expand the 
possible phenomenological studies at, for example, the Large Hadron Collider (LHC) significantly. 

However, the integration of the bremsstrahlung events over phase space hinders further development 
of such generators. Current methods for the automated numerical integration B3]450i of the real matrix 
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elements over the high dimensional phase space require large computer resources. Already for 2 — )• 4 
processes, the obtainable statistical accuracy is a limiting factor for realistic phenomenology. 

On the other hand, shower Monte Carlo programs |[5TH57l tell us efficient algorithms exist to generate 
events closely following the physics 1.58,1 . It is therefore of interest to investigate phase-space generators 
based on especially dipole showers ll59]466l as they employ exact phase-space factorization. A shower 
based NLO phase-space generator uses as a starting point an n-parton final state, which can have additional 
multiple non-partonic particlesQ The origin of this final state does not affect the NLO generator. It can 
be a previously generated Leading Order (LO) event (weighted or unweighted), an event provided by a 
phase-space generator such as Rambo E?!, Sarge [68], Ha AG 1691, Kaleu ITOI or any other source. 
The Forward Branching Phase-Space (FBPS) generator will interpret each of the n parton momenta as a jet 
momentum. Next, following the dipole shower formalism, it performs forward 2 — )• 3 branchings, thereby 
generating the (n + l)-parton phase space. The brancher is constructed in a way that a jet algorithm based 
on a 3 — )• 2 clustering algorithm exactly inverts the branching. As a consequence the observed final state 
is unaltered during the phase-space integration. 

As a result, any observable constructed from the fixed number of jet momenta remains unaltered by 
the FBPS generator and hence does not constrain the real-emission phase-space integration, which can 
thus be viewed as integrating out the partonic degrees of freedom inside the jets. In this sense the jet has 
become opaque and no information can be extracted about the internal jet structure, which is the domain 
of resummed calculations, in particular parton showers. We note that if a dipole shower would use the 
same brancher, the matching of this shower to the NLO calculation becomes a triviality as the shower, just 
as the observables, factorizes from the bremsstrahlung integration. In other words, the shower does not 
alter the NLO event weight nor does it change the partonic jet observable. 

The FBPS generator is an one-particle phase-space integrator, independent from the number of jets. 
Therefore, for a given n-jet configuration, the numerical accuracy is only affected by the number of pos- 
sible 2 — )• 3 branchings, i.e. the number of jet pairs. This means that in order to maintain constant 
statistical accuracy, as the number of jets increase, the Monte Carlo program needs 0{Ncvts x n^) gen- 
erated bremsstrahlung events where A'^cvts is a number of events independent of the jet multiplicity n. An 
additional virtue is that all generated bremsstrahlung events are added to the same virtual event, making 
the infra-red/coUinear cancellations efficient and easy to optimize in the three-dimensional phase space. 
This allows us to use a simple slicing method to facilitate the cancellation of the infra-red/coUinear singu- 
larities. Note that subtraction methods can be trivially implemented as any jet observable does not depend 
on the generated bremsstrahlung event. 

In the first implementation of this method, we use a special 3 — )• 2 clustering jet algorithm, which is 
an augmented 2 — )• 1 jet algorithm. This augmentation adds a recoil parton to the 2 — )• 1 clustering. As 
a result, the NLO jet phase space becomes identical to the LO jet phase space. Specifically, jets resulting 
from the clustering remain massless and the jet algorithm preserves momentum conservation, i.e. particles 
"clustered" with the beam are not discarded. From these observations it is clear, the jet augmentations 
will only modify the last step in the jet algorithm: the clustering prescription. The clustering step is added 
rather ad-hoc to current jet algorithms, which makes it easier to implement the necessary modifications as 
described later on in this paper. 



'At this early stage of developing our method, we do not consider the cases where partonic decays of color-neutral bosons 
may contribute to the n-parton final state. 
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The layout of the paper is as follows. In Section [2] we discuss issues regarding the infra-red safety of 
jet observables. We will propose augmentations of the jet algorithms to make all possible jet observables 
infra-red safe. Because of their enhanced infra-red safe behavior, these augmented jet algorithms have im- 
proved theoretical properties. We will use these properties to construct the FBPS generator in Section [3l 
With this phase-space generator at hand, we build - as a proof-of-principle - a leading-color NLO parton 
level event generator for PP — n gluonic jets. This is detailed in Section IH before we give our con- 
clusions in Section \5\ Finally, two appendices are added. The first appendix details the number and type 
of branchings occurring in the FBPS generator. The second appendix lists the exphcit jet configurations 
utilized in Section |4]to perform the numerical studies. 



2 The Fully Differential Jet Cross Section 

At hadron colliders, jets form the basis of defining the event topology and thereby characterize the un- 
derlying hard scattering event. It is therefore imperative to understand jets in both experiment and theory. 
An essential requirement of the jet defining algorithm is infra-red finiteness, which expresses the fact that 
the addition of arbitrary soft/collinear particles does not alter the final-state jet multiplicity. The infra-red 
finiteness requirement in jet algorithms is by now well understood (see e.g. fTTl '). 

Before the advent of the numerical parton level NLO generators 17211731 . semi-numerical programs 
calculated corrections to differential jet observables. For example, in Ref. (74] the NLO conections to 
the semi-exclusive dijet cross section are calculated for explicitly given values of the dijet mass and the 
rapidities of the two leading jets. This gives a necessarily infra-red safe (finite) correction for each point 
in the dijet phase space. Note that the jet algorithm inevitably forms an integral part of this calculation. 

Current parton level NLO multi-jet generators perform a Monte Carlo integration over the bremsstrah- 
lung phase space independently of any actual jet algorithm. This has the appai^ent advantage that any 
experimental jet algorithm can be numerically accommodated in the Monte Carlo programs. However, 
owing to the infra-red properties of the jet algorithms, the generators only produce infra-red safe results 
for more inclusive jet observables. The resulting predictions from the Monte Carlo generators are not 
infra-red safe for each point in the multi-jet phase space as the LO and NLO jet phase spaces only coincide 
on the boundary defined by soft/collinear emissions. A sufficient amount of jet phase-space averaging is 
required to obtain finite results. For example, one cannot use arbitrarily small bin sizes in representing the 
results of the Monte Carlo integration - sufficiently wide bins are needed. 

As mentioned above while current jet algorithms are infra-red finite, the observables constructed from 
these jets are not necessarily infra-red safe. This is a direct consequence of the clustering procedure 
constituting the last step of the jet algorithms. For example, the dijet azimuthal angle de-correlation is 
a typical non infra-red safe jet observable 1175 117611 . At LO (n = 2), the two jets in the event are exactly 
back-to-back in the azimuthal plane. At NLO, the generation of the jet mass will cause the bremsstrahlung 
events to deviate from the back-to-back configuration leaving uncanceled logarithmic divergences. On the 
contrary, with an infra-red safe jet algorithm, the NLO jet phase space is identical to the LO jet phase 
space and the two jets remain balanced. All that is calculated at NLO is the K-factor; only a third jet will 
induce a de-correlation. 

It is important to note that the Kinoshita-Lee-Nauenberg (KLN) theorem is not a jet phase-space aver- 
aged property. By integrating over all partonic contributions to a. fixed jet configuration, the KLN theorem 
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should already hold. In other words, an infra-red safe jet algorithm must provide a proper cancellation 
between virtual and bremsstrahlung events for each jet phase-space point. As a consequence any jet 
observable constructed through an infra-red safe jet algorithm is finite. 

2.1 Jet algorithms and infra-red safe observables 

To explore the issues with jet algorithms further, we will first look at final-state jets at a lepton collider: 
i'^£~ — Ji • • • J„. Here, a sequential jet algorithm is readily constructed by defining, as a function of 
the cluster momenta {q}, an event resolution measure i?cvt(ci, . . . , Cm). The event resolution separates 
the soft/collinear region from the hard region of phase space. The clusters can be individual hadrons, 
calorimeter cells, tracks and combinations thereof. If i?evt is smaller than the requested jet resolution, 
i?jet, the number of clusters is recombined using the cluster procedure defined by the specific jet algorithm: 
{ci, . . . , Cm} !-)• {ci, . . . , Cm-i}- This is repeated until i?evt > R^et, at which point the remaining clusters 
are identified as the jets. 

Jet algorithms currently used by experiments ( ifTTHTTtiSTIl ') define the event resolution function in 
terms of resolution functions of pairs of clusters: i?evt = minjj Rij : Rij = R{ci, Cj). The minimization 
procedure identifies the least resolved pair of clusters and recombines these two clusters to one new cluster, 
thereby decreasing the total number of clusters in the event by one. In here lies a fundamental issue: either 
the newly formed cluster has a four-momentum, which is massive due to adding the momenta of the two 
clusters Cij = Ci + cj , or overall momentum conservation is violated. From a theoretical point of view, the 
2 — )■ 1 clustering causes the NLO jet phase space to separate from the LO jet phase space. The LO jets 
are massless, while the NLO jets now necessarily are massive. They only match in the exact soft/collinear 
limit, in which case the new cluster is massless: cij = 0. As a result we have infra-red finiteness, but no 
infra-red safety. The necessary care has to be taken when defining jet observables. For a given value of a 
jet observable Cobs, the virtual correction contributes at a single point 6{0 — Cobs) while, because of the 
jet mass, the bremsstrahlung is distributed as log(C' — Cobs)- This behavior is "cured" by allowing, for 
example, sufficient smearing in histogram bins. 

To maintain infra-red safety, we need to both keep massless clusters and maintain overall momen- 
tum conservation when combining clusters. The minimal procedure to do this is by defining the reso- 
lution function in terms of triplets of clusters [82]: i?evt = minjjfe i?(cj, cj, c^). The triplets of clusters 
{cj, Cj,Ck} can be recombined to pairs of clusters {q, Cj}, while maintaining both momentum conserva- 
tion, Ci + Cj = Ci + Cj + Cfc, and keeping the newly formed clusters massless, cf = Cj = oH With this 
type of jet algorithm one can define infra-red safe jet cross sections. As a consequence the fully differ- 
ential cross section dP^^a/dJi ■ ■ ■ dJn, and all possible distributions of jet observables derived from it, 
are infra-red safe. The reason this can be done is that the LO and NLO jet phase spaces exactly match. 
We can therefore construct a phase-space brancher similar to the ones used in dipole showers 116 111641 . By 
choosing the branching map as the inverse of the 3 — )■ 2 clustering used in the jet algorithm, all generated 
bremsstrahlung events are mapped back to the same jet phase-space point. This results in an infra-red safe, 
fully differential jet cross section by virtue of the KLN theorem. That is, both virtual and bremsstrahlung 
corrections contribute to 6{0 — Cobs) only and no smearing is required to obtain a finite result. 

As is clear from the above discussions, it is straightforward to construct a FBPS generator for lepton 
colliders. It would calculate the ii'-factor to a fixed jet phase-space point, i.e. the fully differential jet 

^If flavor-tagged clusters are involved, cf ^ 0, the possible quark mass has to be taken into consideration. 
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cross section. However, for hadron colliders, the incoming partons cause additional complications. The 
current jet algorithms used in hadron colliders augment the lepton collider jet algorithms by including 
a resolution measure of clusters with respect to the beam. If a cluster is combined with the beam, it 
is effectively removed. As a result the remaining clusters violate momentum conservation as we have 
un-clustered momenta. To get infra-red safety, momentum conservation must be preserved during the 
clustering. There are two options: either build up a beam jet, or perform final-state clusterings only. 

The first option is to construct a beam jet: instead of removing the final-state clusters, when combining 
with the incoming beam, they are combined with the respective (separate) beam cluster. Once the event 
resolution passes the jet resolution, we are left with two incoming beam jets and the final-state jets. All 
jets are massless and four-momentum is conserved. However, the two beam-jet momenta are not along 
the incoming (anti-)proton directions. To map onto the LO jet phase space, where the two beam-jet 
momenta are along the incoming (anti-)proton directions, we have to define the jet observables in the 
frame where the two beam jets are along the (anti-)proton directions. That is, we have to perform a 
transverse momentum boost to this frame. From a theoretical point of view, this has the desirable feature 
that the effect of "initial-state radiation" is minimized as this radiation does not affect the observable due to 
the boost. Effectively, the initial-state radiation is integrated out within the jet cone and the KLN theorem 
guarantees a properly defined fully differential jet cross section. 

The second option is to constrain the initial-state clusters to remain along the respective beam direction 
during the clustering phase: the beam particle momenta are only rescaled. While not immediately obvious, 
this can always be accomplished using the 3 — )• 2 clustering mapsH From an experimental point of view, 
this is not a particular desirable option as all radiation is assigned to the final-state jets. 

For a proof-of-principle calculation, the second option is highly desirable as it minimizes theoretical 
complications. It will be used in this paper. As this is a NLO calculation, only one clustering step is 
performed. We start with the partonic scattering paPb Pi - • • Pn+i and reduce this to the jet final 
state JaJb — ^ Ji • • • Jn- Note that a and b are only used to label the incoming partons or jets; no flavor 
information is associated with these labels. Out of the large class of infra-red safe jet algorithms, which 
can be constructed, the explicit jet algorithm used in this paper is as follows: 

1 . Find initial- or final-state parton i and final-state parton j by minimizing the resolution parameter 

Rij = \sij\ = \{±pi+pjf\ (1) 
where "+" is used for i being a final-state particle and "— " for being an initial-state particle. 

2. Given partons i and j of the previous step, find final-state parton k by minimizing 

Rij-k = m.m{Rik,Rjk) ■ (2) 



3. If parton i is a final-state parton: cluster parton i and j, pij = Pi + pj, and use parton k as the 
recoil momentum to make the cluster massless: 

Ji = Pij + (1 - 7)Pfc , 
Jk = IPk 



^This clustering only works for final states with at least one jet. 
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with 7 = 1 + Sij/{sik + Sjk). This maps the three final-state partons onto two massless jets: 
{Pi^Pj-iPk} '-^ {Ji, Jk} while preserving momentum conservation: Ji + J^. = pi + pj +Pk. 

4. If parton i is an initial-state parton, say i = a: cluster the two final-state momenta, pj}^ = pj +Pk, 
and use the initial-state parton a as the recoil momentum to make the cluster massless: 



with 7 = 1 — Sjk/isaj + Sak)- The two final-state partons are now mapped onto one massless jet 
and a rescaled initial-state parton: {pa,Pj,Pk} '-^ {Ja, Jj}< while maintaining overall momentum 
conservation Ja — Jj = Pa — Pj — Pk- 

We will use an inclusive mode of the algorithm. This means, we keep clustering until the desired (LO 
predefined) number of jets is reached. The alternative is to cluster until the jet resolution exceeds the preset 
minimum, after which the event is vetoed, if the number of jets is not equal to the desired number of jets. 
This is the exclusive mode of the algorithm, which is not used in this paper. Note that for reproducing 
the usual NLO n-jet inclusive observables, we have to perform a two-stage run. First, generate the NLO 
if -factors for the exclusive n-jet events, next, add the (n -I- l)-jet events at LO. From this event sample 
the observable can be determined. 

2.2 Defining the fully differential jet cross section 

We want to calculate the fully differential cross section of an n-jet final state characterized by the jet-axis 
momenta Ji, . . . , J„ using the inclusive version of the jet algorithm specified in the previous subsection. 
The jet event kinematics are given by 




'a 



(4) 



XaPa + XbPb = Ji-\ h J, 



(5) 



Jf = 0, 



where denote the incoming hadron momenta and 



•^a,b Pa,b 



(6) 



The collider energy is given by s/S and the momentum fractions Xa and xj, are calculated from the recon- 
structed jets: 



1 



n 




1=1 



(7) 



using the transverse momenta and rapidities of the jets. 




and 




(8) 
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respectively. Note that here we have used the convention p = {E,pz,px,Py). 
We define the differential cross section of the jet observable O as 



da If / \ S^'a 

— = — j d^{xaPa,XbPb^ Jl,---,Jn) dXadXb6[0 -0{Ji,...,Jn)j ^ . . . ^ j ' (9) 

where the jet phase space is given by 



d^{XaPa,XbPb^ Jl,...,Jn) = [j{d^'''^ J, 5{J^) e{Ei)] 5{XaPa + XbPb - Jl Jn 



\i=l 



n '^(^'^^'^ + XbPb-Jl Jn) . (10) 



^1=1 



We consider all jets as indistinguishable and hence have to introduce the "identical-jets" averaging factor 
of (n!)^^. The fully differential jet cross section at LO is given by 



d(-)cTLo (27r) 



4-3n 



dJi---dJn 2XaXbS 

{/a/i,->-/l-/n} 



E FfMFf,{xb) 



M^^\xaPa,XbPb;Jl,---,Jn) 



2 



(11) 



where the Ff^ ^ are the parton density functions of the partons in the beam particles and | A^'''^^^ is the 
squared LO scattering amplitude, spin/color summed (averaged) over final (initial) states. The flavor sum 
runs over all possible, distinguishable partonic subprocesses fafb fi - • • fn that contribute to the jet 
final stated For example, & qq + {n — 2)g final state has n{n — 1) distinct flavor terms to be added for 
one specific initial-state configuration. This way we account for all ways of assigning the distinguishable 
partons of the final state /i, ...,/„ to the jets Ji, . . . , J„. Note that at LO no phase-space integration is 
left for the fully differential jet cross section. However, using our definition, the fully differential jet cross 
section will be symmetric under any exchange of partons without the need of integrating over phase space. 
Once one does the phase-space integration, as in Eq. Q, the usual symmetry factors are recovered. 

Because the jet algorithm preserves explicit momentum conservation and keeps the jets massless, we 
can define a i^-f actor per jet phase-space point Ji, . . . , J^. The NLO corrections to the fully differential 
jet cross section can hence be written as 

dJi---dJn ' ' " dJi---dJn 

In the remainder of the paper we will derive the expression for Knlo and develop, as a proof-of -principle, 
a Monte Carlo integrator for the explicit evaluation of the A'-factor for the pure gluonic contribution of an 
n-jet event at a hadron collider. The K-factor is composed of three contributions, the Born contribution 



''The flavor labels ft denote gluons and massless (anti-)quarks. We omit specifying other than the partonic flavors for reasons 
of keeping the notation simple. For example, we could have a vector boson decaying leptonically in all subprocesses. 
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expressed as "1", the virtual contribution, V, and the bremsstrahlung part, R. We write 

K^l^o{Jl,...,Jn) = 1 + V{jl,...,Jn) + R{Jl,...,Jn) 

= 1 + (^^) . . . , Jn) + . . . , Jn)) , (13) 

where we have factorized the strong-coupling expansion parameter as/2TT and the color factor A'^c- 
Eq. ([T3] ) expresses the cancellation of infra-red singularities per jet phase-space point, V and R on their 
own diverge but their sum gives a finite contribution to the A'-factor. For the calculation of the virtual cor- 
rections, many packages have been developed | l 10 16 22 . 24 - 26] . which can be readily used to calculate 
this part of the K-factoi. On the other hand, the calculation of the bremsstrahlung contribution R requires 
a careful derivation. 

To summarize, for the calculation of a jet observable, we generate, using Eq. the jet configurations 
contributing to the specific value of the observable. For each generated jet phase-space point, we calcu- 
late the LO weight according to Eq. (fTTI ) and the NLO re-weighting multiplicative i^-factor as given in 
Eq. 



3 The Forward Branching Phase-Space Generator 



The explicit construction of the FBPS generator proceeds in several steps. In Section |3TT] the first step is 
taken by the decomposition of the bremsstrahlung phase space into sectors using the event resolution func- 
tion given by the jet algorithm. Each sector is defined through the jet algorithm selecting an unique triplet 
of partons to be clustered. Next, owing to the invertibility of the clustering, we develop in Section [3^ the 
real-emission phase-space formalism based on forward branching off the Born level jet configurations. In 
Section [331 we derive the specific forward branchers, which fill each sector such that the 3 — )■ 2 cluster 
map given by the jet algorithm will recombine the three partons to the same two jets specified by the jet 
phase-space point. Finally, in Section l34l the procedure is validated using the Rambo flat phase-space 
generator. 



3.1 An invertible sector decomposition of phase space 

The bremsstrahlung contribution to the jet cross section, with the jet kinematics specified in Eq. dD, is 
given by 



l-3n 



1 



dJi ■ ■ ■ d Jn I 2 S 



(n + 1)! 



X Ajct{Jl,---,Jn\Pa,Pb,Pl,---,Pn+l) M''°\pa, Pb', Pi, ■ ■ ■ , Pn+l) 



(14) 
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where pa = XaPa and pf, = x^Pb- Compared to the LO case, we now have to sum over all partonic 
subprocesses with one more parton in the final state. As before, of. Eqs. ^ and ([TT]) . the flavor sum 
in combination with the (identical-particle) averaging factor l/(n + 1)! guarantee the correct symmetry 
factors for the (n + 1) -particle final states. 

The generalized jet delta-function Ajct decides whether a specific subprocess with its parton kinematics 
contributes to the bremsstrahlung factor R at the jet phase-space point Ji, . . . , J„. This jet delta-function 
is equal to unity if the jet algorithm clusters the parton momenta pa,Pb,Pi, ■ ■ ■ , Pn+i to the jet-axis four- 
vectors Ji, . . . , J„ and is zero otherwise. The function (as used here) is flavor blind - the jets as well 
as the partons are indistinguishable, therefore, Ajet has to be symmetric under any exchange of jet and 
parton momenta. To integrate over the jet delta-function, we make use of the 3 — 2 clustering algorithm 
discussed in the previous section. This allows us to expand the jet delta-function over a sum of dipoles, 
each selecting three partons, which will be clustered by the jet algorithm to two jets and, as a result, a 
value for the resolution parameter Rij-k will be returned for any of these combinations. The expansion 
can then be written as follows: 

Ajet(-^l, • ■ ■ ,Jn\ Pa,Pb,Pl, ■ ■ ■ ,Pn+l) 
= 6{Xa - Xa) 6{Xb - Xb) X 

Qveto 

^3et{Ji;Jj\pkvPk„+i;Pkj) — ' ' n ^^J^-Pks) 

iT^j fcG5„+i s^i,j 
i,j=l,...,n s=l,...,n 



+ 6{xb-Xb) Y ^}<^^{'Ji\Pci'^Pkr.+i^Pk?) Q7k°+rM n '^(•^^ 
fce5„+i s^i 
i=l,...,n s=l n 



Pks 



+ 6{Xa-Xa) Y ^3^ii-^i\Pb''Pku+i^Pk,) QlkT+iM H ^i-^s - Pks) , (15) 
fceS„+i s=ii 
i=l,...,n s=l,...,n 

where the sums are over all possible pairs of jet-axis momenta Ji , . . . , J„ and permutations of the brems- 
strahlung four-momenta pi, . . . ,Pn+i- The index vector k describes the elements of the permutations 
Sn+i of the set {1, 2, . . . , n + 1}. The permutation sum ensures that all dipole-spectator configurations 
and their respective (n — 2)! phase-space combinations in the leftover parton momenta are taken into 
account. 

The sector veto (or, just as well, jet resolution) cut QJ^^j implements the first two steps of our jet 
algorithm proposed in Section IZTl According to Eq. ([T]), we evaluate 

-Rmin = Rvw = min (16) 

vw 
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where v^w = 1, . . . ,n + 1 and v < w. If we take v,w^r = l,...,n + l and define 

i?|-nin = min Rfju,.r = miii mm{Ri)r , Rwr) , (17) 

r r 

cf. Eq. Q, we can formulate this cut as 

®im;j — ^ {Rm'm Rim) ^ (^min Ri^\j) i (18) 

using the Heaviside step function 6{x), which equals one if a: > and is zero otherwisej^] Step 3 and 4 
of the proposed jet algorithm are executed when one computes whether the generic Ajct building blocks 
are different from zero. They will return one, only if the parton momenta reconstruct to the two given 
jet momenta or, in the initial-final state cases, to the one given jet momentum. More specifically, we 
understand Ajct(J; J' \ Pi,Pm',Pj) = 1, if and only if J = pim + (1 — "y)Pj and J' = jpj, cf. Eq. 
Neither in the vice versa case J -H- J', nor in any other combination we find the generic Ajct / 0. This 
way we avoid double counting when we permute over all parton momentum configurations. The emitter- 
emitted parton symmetry in 0™J.° , however, leads to double counting the same event in the final-final 
parton sum of Eq. (fT5l) . which we remove by multiplying the factor 1/2 to the vetoH Since the order of the 
bremsstrahlung momenta is permuted, we guarantee that all combinations are tested (in Ajct as well as the 
product of the 6{Js — PkJ terms) to make sure that the bremsstrahlung events are selected, which match 
the considered jet phase-space point kinematics and, therefore, give a contribution to R{Ji, . . . , Jn). 

We observe in Eq. (fT5l) that the jet delta-function breaks up phase space in two types of sectors: final- 
final state sectors and initial-final state sectors. Note that in principle there could be an initial-initial state 
sector as well. However, this will only occur if we allow for the build-up of beam jets. 

We finally note that for each jet phase-space point exactly one sector contributes. With R'^^^ (based 
on i?min) a global event resolution measure is given, which depends on all initial-state and jet momenta. 
This partitioning of phase space is dictated by the event resolution function given by the jet algorithm, see 
Eq. ([T])- As a consequence the phase space is invertible: given the jet four-momenta Jj and Jj one can - 
by inverting the cluster map of the jet algorithm - generate the three-parton configurations for each sector, 
which will cluster back to these two initial jets. Such a forward branching Monte Carlo integrator exactly 
integrates out the internal jet structure. In the next section we will formulate these forward branchers. 



3.2 Phase-space construction through forward branching 

To construct the forward brancher for a sector, we have to integrate Eq. ([T4l ) over the jet delta-function of 
Eq. (fT5l) . The jet delta-function selects those (n + l)-parton final states, which reconstruct to the given 
n-jet phase-space point. We can turn the approach around and use the jet delta-function as a prescription to 
explicitly generate the n+ 1 bremsstrahlung parton momenta given the n-jet momenta. This establishes the 
forward-branching picture, which in addition allows for the avoidance of the dipole and permutation sums 
of Eq. (fT5] ). To see this, we can write down the final-final state piece of Eq. (fT4l) for a single subprocess 



^The jet clustering of Section 12.11 is symmetric under the exchange of the emitter and emitted parton (Rwv ~ Rvw and 
Rwv;r ~ Rvw;r)- Therefore, we do not have to consider to > ?; to determine the iimin and R'^^i^- 

^Later on, we resolve this issue by partitioning the phase space further according to the different parton emitter settings. 
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neglecting all prefactors, including the Bom matrix element: 

Rpf{Ji, • • • 1 Jn) 

Yl Yl d^{Jij^Ph, Pk„+i , Pk, ) Ajet {J^ ; Jj I Pk, , ; Pk, 

i+j A:G5„+i 
i,i=l,...,n 



X 



Qveto 



(19) 



where Jij = Ji + Jj . We observe that the phase-space integration has been broken up into many factorized 
pieces of splitting dipoles. Furthermore, we can exchange the order in performing the dipole phase-space 
integrations and summations. The integration over the three-parameter phase spaces can be accomplished 
through Monte Carlo techniques. We can treat the explicit dipole and permutation sums similarly: in- 
stead of carrying them out, we can choose dipole and parton configurations at random^ We just have to 
keep track of and include possible weights that may occur in the selection of dipoles and bremsstrahlung 
partons. 

Eqs. (fTTl) and (fT4l ) have complete flavor sums running over all possible subprocesses that contribute to 
the n-jet and (n + l)-jet final states, respectively. We can maintain this inclusive-flavor approach in the 
forward generation of the real-emission events. No knowledge of the particular LO process and its flavors 
is needed apart from the given set of the jet-axis momenta, which we interpret as the initial four- vectors 
before the parton branching. The forward branching occurs, in principle, independently of flavor; the 
pure generation of the bremsstrahlung momenta in fact has no flavor dependence. The only place where 
flavor conditions enter is in combining phase space with the matrix element for the randomly chosen 
subprocess containing 2 — )• n + 1 strongly interacting particles: the number of 3 — )• 2 clusterings as given 
by combinatorics may reduce owing to flavor constraints If] 

To simplify the discussion, we focus on the pure gluonic case. Consequently, the flavor sums in 
Eqs. (fTTl) and (fT4l) collapse to single terms. We also can simply arrange to set the first n partonic four- 
vectors according to the jet-axis momenta. In the final-final case for example, two of these, the emitter 
and spectator momenta, will change owing to the generation of the additionally emitted parton, which 
we can always choose to label by n + 1. For ordered amplitudes, one may insert the new parton right 
after the emitter parton / — 1 and shift all subsequent ones by one, / — )• / + 1. Because of the forward 
construction of the parton momenta, the criteria underlying the generic A^ct {Ji', J j \ Pi,Pr]Pj) terms will 
be satisfied by construction. Thus, these terms are redundant and the constrained generation of pi, pr and 
Pj already accounts for the (n + 1) n (n — 1) combinations of arranging three partons to be clustered to the 
two jets picked for forward branching. Still, with respect to the non-branching part of the final state, we 
have (n — 2)! possibilities to assign the parton momenta with certain jet momenta. However, owing to the 



^Owing to the flavor blindness of the jet definition, the parton configurations have to be varied too as long as the chosen 
subprocess contains different parton flavors. We can use any set of indistinguishable particles though to reduce the initial number 
of (n + 1)! possibilities. 

^As an alternative the FBPS generator may be designed such that the parton flavors are treated as in dipole showers. For 
example, select a flavor assignment at LO as in Eq. i ll It . now consider an initial-state branching; if the LO subprocess has an 
incoming quark, there are two bremsstrahlung contributions: q ^ qg and g ^ qq where the gluon and anti-quark are radiated 
off, respectively. The corresponding n + 1 matrix element then determines the weight of the selected option. 
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symmetry of the final state, the leftover permutation sum can be replaced by a multiplicative factor. The 
QJ^^j is included to the three-parton phase space and acts as a phase-space cut implementing the jet res- 
olution criteria. Applying similar arguments to the initial-final state cases, we can take all modifications 
and write the result of the phase-space integration as 



i?(Jl, . . . , Jn) 



1 



1 



(27r)3 n(n + l) 



(Ja = XaPa, Jb = XbPb] Jl, 



) Jn) 



E 

ij=l,...,n 



n — 1 



d^^cto{Ji,Jj ^Pi,Pr,Pj) M^^\ja,Jh;{Ps = Js}s^ij,Pi,Pr,Pj) 



+ ^ d^yeto{Ja,Jj'-^Pa,Pr,Pj) 



XaFfJXa) 



M^^\pa,Jb;{Ps = Js}s^j,Pr,Pj) 



+ ^ d^yetoiJb,Jj ^Pb,Pr,Pj) ^ p . N 



M^'^\ja,Pb;{Ps = Js}s^j,Pr,Pj) 



(20) 



Note that the (n+ 1) ! term in the denominator has been combined with the multiplicative numerator factors 
(n — 2)! and (n — 1)! for final-final and initial-final state branchings, respectively. Also, r = n + 1 and 
Pa,b = Xa,bPa,b whilc b = Xa,bPa,b- As bcforc, the dipole sum and the three-parameter phase space can 
be calculated using a Monte Carlo integration, i.e. / — > (I/A'mc) Yl,i=i Nmc d P^^^ ■ 

The dipole factorization of phase space is obvious from the equation above. The 2 — )• 3 differential 
phase-space volumes can be described by dipole or antenna phase-space factors, which are also used in 
shower algorithms. In our calculation the combination with the matrix element then fully specifies the 
branching, which in showers is achieved only approximately by the use of the splitting function. The 
final-final state antenna phase space, cf. ||72][83], is given by 



d^ycto{Ji,Jj ^Pi,Pr,Pj) 



©veto -I 
ir-j TT L 



r) r) ^ ^i'l' ^ ''^fj r) 



(21) 



whereas the initial-final state antenna phase space, cf. II731I84I is expressed as 

P ■ J 



d^wetoi^Ja — XaPai Jj ' ^ XaPaiPriPj) — ®ar;j 



1 d^^^pr 



(22) 



27r 2Er \Pa-Jj-Pa-Pr 

In the latter case, pj is given by momentum conservation, i.e. XaPa — Jj = XaPa — Pr — Pj- The new 
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momentum fraction Xa can be calculated using the condition p'j 



= 0; we obtain 



Jj ' Pr 



(23) 




In both cases above, the phase-space factors have to be supplemented by the corresponding sector veto 
(or jet resolution) cut. The cuts guarantee the partitioning of the bremsstrahlung phase space such that no 
overlapping regions can emerge. The actual sector phase-space volume is then measured by the Monte 
Carlo integration by means of the sector veto cut. 

The formulation of the effective FBPS generator used in the remainder of the paper is now simple. 
Given an n-jet event, we generate the bremsstrahlung events in the following manner: with probability 
(n — l)/(n + 1), select a pair of final-state jets randomly and perform a final-final state branching. 
Else, select one of the two incoming partons and one final-state jet and perform an initial-final state 
branching. Appendix lAl gives more explanations regarding this selection. The bremsstrahlung event now 
has n + 1 final-state partons, which reconstruct back to the original jet configuration using our specific 
jet algorithm. We repeat the procedure until a sufficient number of bremsstrahlung events have been 
generated to estimate the ET-factor for this particular jet event. This in fact is the execution of the Monte 
Carlo integration, whose uncertainty can be controlled by the number of generated Monte Carlo events 
per jet phase-space point. 

3.3 Forward branchers 

To completely assemble the forward-branching buildup of the bremsstrahlung phase space, we still have 
to define the generated three-parton final state in terms of the original jets and the dipole phase-space 
integration variables. In doing so we have to respect the constraint that the jet algorithm clusters the 
three generated partons back to the two initiator jets. In the next two subsections, we will formulate the 
branchers, explicitly designed for this task. 

3.3.1 The final-final state brancher 

From the final-final state phase-space factor of Eq. (|2T]) we extract the phase-space factor for the FBPS 
generator. We write 



where the sector veto cut has been introduced in Eq. ([TSl l. Several comments are in order. We have defined 
Uki = Ski/sirj with Ski = {Pk + Pi)^ and Sirj = {pi +Pr + Pj)"^ = {Ji + Jj)"^- Apart from the kinematic 
constraint, yij. + yrj < 1, we added an additional constraint, y^j < yij, which divides each sector into two 
(sub)sectors, breaking the emitter-emitted parton (i o r) symmetry. As a result, the factor 1/2 formerly 
present in Eq. (|2TI ) is dropped here. The additional constraint is needed to accommodate for the integration 
over asymmetric functions in pi and pr, e.g. over ordered amplitudes (or over quark-gluon states, which 
is important for later applications). The distinction is particularly important when we combine the sectors 
with ordered matrix elements as each sector has its own singularity structure. It manifests the notion of 



d^vcto{Ji,Jj ^Pi,Pr,Pj) 



- Sirj dyir dyrj — e{l-yir- yrj) {yij - yrj) &Jr;j (24) 
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Pi, pr and pj respectively being the emitter, emitted parton and spectator in the phase-space branching. 

The algorithmic description of the final-final state brancher is as follows: starting from two final-state 
jets Ji and Jj; 

1. generate the integration variables yir, Vrj and (j)/ {2tt) on the interval (0, 1] and fulfill the constraints 
Vir + Vrj < 1 and t/rj < i/ij. Noticc that the sector veto cut 0™.*° guarantees y^r < Urj, cf. Eq. ([TSl l. 

2. rescale the four-momenta Ji and J, : 



J ki = Ji + jJj, 
I % = (1-7) Jj , 
where 7 = y^r such that we find kf = yir Sirj. 



(25) 



3. determine the final-state partons pi and pr by invoking the phase-space decay of ki ^ pi -\- pr with 
the on-shell condition pf = Pr = 0. 

4. set Pj = kj. 

5. pass the event if and only if the sector decomposition cut 0™*/ has been satisfied. Assign the weight 
Sirj 7r/2 to the eventj^ 

Using this construction procedure, the final-final state clustering of the jet algorithm maps the partonic 
set {pi,Pr,Pj} created through this FBPS generator back onto the jet pair { Jj, Jj}. Note that each of the 
forward branchers would generate the "same" one-particle phase space, if the sector veto cut G™.*? was 
removed from Eq. (l24l ). An additional integration over the n-jet phase space would generate the whole 
{n + 1) -particle bremsstrahlung phase space. 



3.3.2 The initial-final state branclier 

From the initial-final state phase-space factor of Eq. ((22]) we extract the corresponding phase-space factor 
for the FBPS generator: 

Because the initial-state parton is distinct from the final-state partons, no additional ordering requirements 
are present. Note that - as in the previous case - by removing the sector veto cut 6^$:*^, one allows for 
this forward brancher to generate the whole (n + 1) -parton bremsstrahlung phase space. The algorithmic 
construction underlying this FBPS generator is outlined below. It is set up such that the initial-final state 

'This may be supplemented by a possible weight from the generation of the integration variables. For example, one may 
rewrite da as s ci(log s), which would generate an additional weight to be included. 
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clustering of the jet algorithm maps the generated triplet of momenta {xaPa,Pr,Pj} back onto the initial- 
final jet pair {xaPa, Jj}- The generation of the momenta then proceeds as follows: starting from an 
initial-state parton XaPa and a final-state jet Jj; 

1 . generate the one-particle phase-space momentum pr within the appropriate integration boundaries. 

2. having generated pr, calculate 

Z = Xa/Xa = 1 + (Jj ■ Pr)/[xaPa ■ (Jj - Pr)] , 

(27) 

Pj = -Pr + {Z-1) XaPa ■ 

3. pass the event with weight (^p-pjr^^ assigned, if and only if the jet resolution cut @ar^.j has 
been satisfied. 



3.4 Numerical validation of the phase space generator 



We want to verify the FBPS generator on itself, before we use the generator to calculate the i^-factor 
for an n-jet phase-space point. For this purpose we do not add in the contributions stemming from the 
matrix elements and PDFs. This is an important validation to ensure the correct treatment of the weight 
generation during the build-up of the bremsstrahlung phase space. We use the flat phase-space generator 
Rambo ||67l for the numerical validation of the FBPS generator. The phase space generated by Rambo, 
dTZn, is connected to the customary flat phase space d ^ through: 



d^(XaPa,XbPb l-> Jl, . . . , Jn) 




XaXbSY^-^^ 

-l)!(n-2)! 



dnn{XaPa,XbPb H> Ji, . . . , Jn) . 



(28) 



For this test, we do not include parton density functions; instead we choose the parton fractions uniformly 
between zero and one. We define an n-jet event at collider energy of 7 TeV using the following selection 
criteria for the n jets regarding their transverse momentum, rapidity and geometrical jet-jet separation: 
Pj) > 250 GeV, \yi\ < 2.0 and ARij > 0.5. We present the results of our tests in Figures [T] and [2] where 
we exemplify the FBPS validation by means of comparisons of distributions for two distinct observables. 
Here, we define 



Ht = \ 



E 



-(0 
Pt 



s 



(29) 



as a dimensionless (scaled) variant of an ordinary Ht variable where instead of using the scalar sum of 
the jet pts, the squared quantities have been summed up. With the second observable, we want to look 
into an angular distribution, namely the azimuthal angle between the two leading jets. This variable is 
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Figure 1: The Ht and A(f>i2 pure phase-space comparisons between clustered {n + l)-particle Rambo events 
(blue lines) and forward-branched n-particle RAMBO events (red lines) for n = 2-jet (upper graphs) and n = 3-jet 
(lower graphs) configurations. The lower panels show the ratio between the clustered and branched predictions 
minus one. The black lines represent the results of the corresponding n-particle Rambo generations. To focus on 
the phase-space validation, the matrix-element and PDF weights have not been included in these calculations. 
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calculated by 



A(pi2 = arccos 



( 




) 



(30) 



For various n-jet final states, the Ht and A(j)i2 distributions resulting from our phase-space tests are 
displayed in Figures [T] and [2] All plots contain three curves labeled by "Rambo", "Rambo+Branching" 
and "Rambo+Clustering" (represented by the black, red and blue lines, respectively); in the corresponding 
lower panes, we always show the ratio, subtracted by one, between the latter two predictions. In general 
we observe steeper tails in the Ht spectra for an increasing number of jets. The A(j)i2 distributions show 
the opposite behavior; their low angle-difference bins become more populated owing to a larger amount 
of phase space being filled. For the same reason of enhanced phase-space filling, the total energy is shared 
among more jets, which leads to the suppression of the Ht tails. 

We now explain the three different predictions, which we use for the validation and show in the figures. 
The "Rambo" curves are obtained from jet momenta generated according to Eq. (1281 ) where we define the 
"LO" n-jet phase space as the n-particle uniform phase space. The jet momenta satisfy the acceptance 
cuts given above. To produce the "NLO" n-jet phase space, we generate n + 1 particles in flat phase 
space with the help of RAMBO and apply our jet algorithm to find n jets from which we can calculate 
the jet observables. Again, these jets have to fulfill the acceptance criteria. This procedure gives us 
the "Rambo-i-Clustering" predictions in Figures [T] and |2] The differences seen between "Rambo" and 
"Rambo-i-Clustering" visualize the pure phase-space effect when generating the "NLO" corrections. 

We can now validate the FBPS generator: we construct "LO" n-jet configurations using the flat phase- 
space generator as we did for the "Rambo" predictions. For each configuration, we subsequently generate 
bremsstrahlung events using the FBPS. As the generated events always reconstruct back to the origi- 
nating n-jet event, we only have to average over the generated event weig This determines the 
"Rambo-i-Branching" curves, which have to coincide - apart from statistical fluctuations - with the re- 
spective curves of the "Rambo-i-Clustering" procedure. The ratio plots in Figures [Hand |2] illustrate how 
well this is achieved by directly comparing the flat phase-space "NLO" predictions and the FBPS gen- 
erated "NLO" predictions. We see excellent agreement between the two results and thereby validate the 
FBPS generator. It is interesting to note that at "NLO" the A4>i2 distribution for two jets does not show 
the usual feature of de-correlating. Because the jet algorithm is infra-red safe for all jet observables, the 
two jets are always exactly back-to-back (in the azimuthal angle) for both "LO" and "NLO", i.e. this dijet 
observable is not affected by the initial-state radiation. 

4 The Gluonic Jet Generator 

As a proof-of-principle we describe in this section a generator, which calculates the NLO K-factor per 
jet phase-space point for the pure gluonic part of n-jet production at hadron colliders in the leading-color 
approximation. 

The determination of the ii'-factor for a given n-jet event requires the evaluation of a single virtual event 
and a three-dimensional Monte Carlo integration over the bremsstrahlung phase-space sectors defined by 

'"Of course, we also cross-checked that the backward clustering indeed recovered the n-jet phase-space point. 
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Figure 2: The Ht and A(f)i2 pure phase-space comparisons between clustered {n + l)-particle Rambo events 
(blue lines) and forward-branched n-particle Rambo events (red lines) for n — 6-jet (upper graphs) and n = 10-jet 
(lower graphs) configurations. The lower panels show the ratio between the clustered and branched predictions 
minus one. The black lines represent the results of the corresponding n-particle Rambo generations. To focus on 
the phase-space validation, the matrix-element and PDF weights have not been included in these calculations. 
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the jet resolution scale. This is discussed in Section |4~T] Finally, in Section l4!2l our results are presented: 
we verify the cancellation of the slicing parameter in the K-factor and the evaluation of the full leading- 
color NLO corrections up to 15-jet configurations. 



4.1 Evaluation of the X-factor 



In terms of ordered amplitudes m, the fully differential gluonic leading-order, leading-color contribution 
to the n-jet cross section is given as 



(2vr) 



4-3n 



d Ji ■ ■ ■ d Jn 



2 XaXf) S 



E 

(T e 5,1+1 



Fg{xa)F,{xi,)N^{N^-l) 



fffWfj- P J J J ' 



+ o 



(31) 



where the color, A'^c, dependence has been made explicit. The sum is over all (n + 1)! different permuta- 
tions of the ordered amplitudes and the jm^*^)!^ are the corresponding helicity averaged/summed squared 
matrix elements. Note that the take values from 6, 1, . . . , n where = x^P}). We can now define the 
fully differential NLO cross section through ordered /c-factors. 



Cnlo 



(2vr) 



4-3n 



IJl - ■ ■ dJn 



2 XaXb S 



Fg{xa)Fg{x,)N^{Nl-l) 



^ ^ ^NLO {^Jai Ja\i Ja2T ■ ■ ■ ^ JcTn+i ) ^ (^Jai JaiT Ja2T ■ ■ ■ ^ Ja^+i ) 



2 / 1 



(32) 



with Ja = XaPa- Noticc that each ordering has been assigned its own /c-factor. 

We further detail the ordered A;-factor similar to Eq. ([T3] ) by dividing the virtual contribution into two 
parts, which we call Vo and /. The former is proportional to the LO term and contains the singularities; 
the latter describes the finite virtual corrections. We now furthermore specify helicity dependent fe-f actors 
and write representative for all orderings 



Vai^Ja, Jb, Jl, ■ ■ ■ T Jn) + f{^JaiJbiJli---^Jn) + f i^Ja'^ i ^ i ■ ■ ■ } Ji 

where 

f( T'^a T^^ 7^^"'\ 

J ' "^fe ' "^1 1 ■ ■ ■ 1 -^n j 



f( 1^'^ T-^^ T'^"^ 



m(o)(J^,J,^^J,^^...,J^) 



(33) 



(34) 
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Note that the Vo part is independent of the helicity structure and is known in analytical form II72[|73II . The 
helicity information is carried by the LO amplitude. We use the code of Ref. (5) to calculate the finite, 
helicity dependent contribution /. The unresolved soft/collinear contributions contained in r are obtained 
using a phase-space slicing method II721I73I . Similarly to the singular virtual terms, this analytically cal- 
culated contribution is helicity independent and proportional to the Born level expression. We add it to 
the Vjy term, which now becomes dependent on the slicing parameter: Vo — )■ foCsmin). while the resolved 
contribution is contained in r — r(smin)- The ordered resolved bremsstrahlung factor r(sniin) is defined 
such that it is independent of the helicity choice. This is done in the same manner as for the unordered 
bremsstrahlung factor R = in Eq. (|20l ). We define the bremsstrahlung factor as the ratio of the 

helicity summed bremsstrahlung contribution divided by the helicity summed LO contribution, schemat- 
ically written as |m^,Vmf"|2 ^ ( E jKrems""^' I V kS'""""!') x |m,^o''-'^'f _ By doing so we 
have arranged for having all the helicity dependence carried by the jets Ji, . . . , J„, i.e. the Born level 
ordered amplitude. This approach of defining the hehcity dependence gives us the correct behavior: the 
soft/collinear limit is found to be proportional to the LO helicity dependent ordered amplitude. Further- 
more, if the helicity sum over the jets is carried out, one retrieves the helicity summed bremsstrahlung 
amplitude. 

At this level of the Monte Carlo program development, it is convenient to have an explicit log^(smin) 
dependence in the unresolved parton contribution. In the previous section we validated the FBPS generator 
based on and using flat phase-space generation. However, the soft/collinear limit is hardly probed by 
distributing the momenta uniformly in phase space. As the Smin dependence of the Vq part has to cancel 
against that of the FBPS generator, we get an excellent probe on the crucial correctness of the soft/coUinear 
behavior of the FBPS generator. 

In the future we can switch to a subtraction method, eliminating the dependence on the slicing param- 
eter. For the class of infra-red safe jet algorithms, this is almost trivial. The observable jet final state 
is invariant under the bremsstrahlung Monte Carlo integration. This means we only have to add to the 
/c-factor, an integral over the unresolved phase space. The integrand is simply given by the difference 
between the bremsstrahlung matrix element squared and its antenna approximation. The Smin parameter 
then becomes equivalent to the so-called a parameter introduced in Ref. 1851 . 



4.2 Numerical studies of NLO high multiplicity jet events 

For all numerical studies we use single, exclusive n-jet events. As for the numerical validation of the 
FBPS generator in Section [331 all events pass the jet cuts p;^ > 250 GeV, \yi\ < 2.0 and ARij > 0.5 
at a collider energy of 7 TeV using the CTEQ6M PDF set |[86l . The renormalization/factorization scale is 
set to one half times the average dijet mass. This scale choice is closely connected to shower Monte Carlo 
approaches where the dijet mass is often used as the starting scale for branchings off the particular dipole 
antenna. 

We first have to study the dependence on the slicing parameter Smin for explicit jet configurations. To 
this end we calculate the resolved, r(sinin)> and the unresolved, Vj){snim), helicity independent contribu- 
tions to the /c-factor. The results are shown in Figures [3l |4] and [5] together with the helicity independent 

part of the /c-factor, 1 + ^ ° ) [''(■^min) + ^^D(smin)], for a single 2-jet, 3-jet, 4-jet, 5-jet and 8-jet 
event and two different orderings per kinematic configuration. The notation in the figures is such that 
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Figure 3: The cancellation of the Smin dependence for the 2-jet and 3-jet configurations given in AppendixiBl The 
upper panels show the resolved contribution in red, the unresolved contribution in black and the helicity independent 
part of the fc-factor defined in Eq. (l33T l in blue. In the lower panel, the dependence on s,„in of the helicity independent 
part of the fc-factor is shown in finer detail. The different color orderings (indicated in the respective plot titles) use 
the same events in the Monte Carlo evaluation, leading to correlated statistical fluctuations. 
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Figure 4: The cancellation of the Smin dependence for the 4-jet and 5-jet configurations given in AppendixiBl The 
upper panels show the resolved contribution in red, the unresolved contribution in black and the helicity independent 
part of the fc-factor defined in Eq. (l33T l in blue. In the lower panel, the dependence on s,„in of the helicity independent 
part of the fc-factor is shown in finer detail. The different color orderings (indicated in the respective plot titles) use 
the same events in the Monte Carlo evaluation, leading to correlated statistical fluctuations. 
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0.00552438 


1.09± 0.05 
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243±2 
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1.18±0.08 
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1.10± 0.08 
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392±3 


125.234 


1.30± 0.13 


0.282847 


1.17± 0.13 


0.0498892 


1.18±0.13 


5 


366±4 


5941.55 


0.94± 0.17 


849.054 


0.87± 0.17 


31.5083 


0.80± 0.17 


6 


529±5 


1202.54 


1.15± 0.24 


69.0066 


1.06± 0.24 


0.469815 


0.82± 0.24 


8 


650±7 


26732.0 


1.41±0.34 


1364.49 


1.32± 0.34 


1.41604 


1.15± 0.34 


10 


844± 11 


6575.23 


1.49± 0.49 


579.066 


1.26± 0.49 


6.09232x10-6 


0.97± 0.49 


15 


1264±20 


4690.02 


1.39± 0.95 


671.554 


1.28± 0.95 


4.37178x10-^ 


1.24± 0.95 



Table 1: The LO ordered amplitude squared jm^^-' ( Ja, Jf,, Ji, . . . , J„)p and its corresponding r(sinin) and ordered 
fc-factor as defined in Eq. (|33] | for an exclusive n-jet event. The explicit jet momenta for the different jet multiplicities 
are given in Appendix iBl The slicing scale Smin is set to lO^'' x S and the Monte Carlo integration over the 
bremsstrahlung phase space has been done with 100,000 generated events. 



m^^\l, 2, . . . , n + 2) = m^'^\Ja, Jb, Ji, . . . , Jn) with l-H-a, 2-H'6, 3^1 etc. For each antenna, we 
only have to perform a three-parameter integral, giving us good control over the cancellations. The graphs 
demonstrate that the cancellation of the Smin parameter dependence is achieved in a satisfactory manner at 
values of the order of 10^^ x S and smaller. Moreover, we were able to maintain good numerical stability 
down to values of Smin < 10^^ x S {= S'coUider in the figures), even though we did not use any adaptive 
Monte Carlo integration such as VEGAS [87] to obtain these results!^ 

We can now proceed to make some predictions for the n-jet configurations listed in Appendix |B] In 
Table [T] we show the results for several helicity configurations of the calculations of the /c-factors de- 
fined in Eq. (|33] ) multiplying the ordered amplitudes squared. As can be seen from the table, the scale 
choice described above leads to relatively small fc-factors, i.e. the normalization of the LO prediction is 
fairly close to the NLO rate. In other words, the correction of the LO weight due to radiative correc- 
tions is of the order of one. In addition to the A;-factors, our table also displays the numerical results 
''(smin) of the brcmsstrahlung phase-space integration using the FBPS generator. As expected, for a fixed 
number of bremsstrahlung events, the uncertainty on the integration results scales with the square root 
of the number of dipoles, i.e. it grows linearly in the number of final-state jets (except for the two high- 
est jet multiplicities, where with the chosen number of events we did not achieve a sufficient accuracy 
for determining a reliable uncertainty estimate). The number of final-final and initial-final dipoles is 
n(n — 1) and 2n, respectively. To obtain an integration uncertainty independent of the n-jet multiplicity, 
we would have to scale the number of bremsstrahlung events evaluated in the Monte Carlo integration as 
n(n — 1) + 2n = n(n + 1). 

From Table \T\ one also reads off that the relative errors on the /c-factors become sizeable for a large 
number of jets. This is when the values of the r-factors turn big and one finds large cancellations between 



"For future applications, the resulting three-parameter integration can readily be optimized by important sampling and adap- 
tive stratification. 
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Figure 5: The cancellation of the Smin dependence for the 8-jet configurations given in Appendix |B] The upper 
panels show the resolved contribution in red, the unresolved contribution in black and the helicity independent part 
of the fc-factor defined in Eq. ( |33] ) in blue. In the lower panel, the dependence on Smin of the helicity independent 
part of the fc-factor is shown in finer detail. The different color orderings (indicated in the respective plot titles) use 
the same events in the Monte Carlo evaluation, leading to correlated statistical fluctuations. 



the single terms in Eq. (1331 ). Owing to the use of a slicing method and its explicit dependence on log(smin)^ 
we are not able to avoid this behavior easily without switching to an approach based on subtractions. 
Because of these large cancellations, the absolute values of Ar multiplied by ^^f^^ ~ 1/20 are hence 
pivotal in determining the uncertainty on the /c-factors. 



5 Conclusions 

In this paper we derived a new type of NLO phase-space generator. This forward-branching phase-space 
generator has the property of inverting the clustering occurring in the jet algorithm. Because of this, the 
bremsstrahlung phase space of a fully exclusive jet final state is generated. The bremsstrahlung events 
constructed in this way do not change the value of a jet observable; they are all added to the single virtual 
contribution. This gives a perfect cancellation of the divergent pieces - the soft/collinear real and virtual 
contributions - as dictated by the KLN theorem for any jet observable. 

However, the current jet algorithms used by the experiments employ a 2 — )• 1 clustering scheme. 
Furthermore, a beam jet is not defined leading to transverse momentum imbalance in jet events. As a 
result the LO and NLO jet phase spaces are different and only match at the soft/colhnear boundary. This 
makes the jet observables infra-red finite, but not necessarily infra-red safe. The KLN theorem becomes 
applicable only after some phase-space averaging. As a consequence, the fully exclusive multi-jet event 
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is not defined at NLO. 

Any infra-red finite jet algorithm using a 3 — )• 2 clustering scheme, which includes beam jets, suffices 
to define the fully exclusive multi-jet differential cross section. For the purposes of this paper, we worked 
out a theoretically simple jet algorithm such that fully exclusive jet final states can be defined at NLO. 
Nevertheless, the difference between an observable determined by this more theoretical jet algorithm and 
one currently used by the experiments amounts to a finite correction. This correction is readily calculated 
for a particular observable by an extra bremsstrahlung phase-space integration, extending the applicability 
of the forward-branching phase-space generator beyond jet algorithms using 3 — 2 clustering schemes. 

Since the forward-branching bremsstrahlung phase-space generator does not alter the jet configuration, 
a K-factor can be defined for a given multi-jet phase-space point. The i^-factor can be determined from 
the leading-order probability associated with the particular jet phase-space point. In this sense, the proba- 
bilistic interpretation of the NLO prediction is restored by quantifying it as a positive-weight adjustment. 
We validated our method in two steps: first, we verified the construction of the bremsstrahlung phase 
space; second, we calculated the radiative corrections for individual events with up to 15 jets. As our test 
scenario, we chose 2 — )• n gluon production in the leading-color approximation. 
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A The Enumeration of Branchings 

We can choose the branching pair of jets in Eq. (l20l) by Monte Carlo means. However, it is useful to 
control the ratio of initial-final and final-final branchings for optimization purposes. This ratio is given 
by the number of ways the n + 3 partons can be clustered to n jets plus 2 incoming partons. If no flavor 
constraints are taken into account, the derivation of the ratio for the jet algorithm used in this paper goes 
as follows: 

- Given an initial-state parton, there are n (n + 1) possible clusterings since it is the final-state parton 
pair that effectively gets combined. This number is based on treating all final-state partons as dis- 
tinguishable particles. The recoil is taken by one of the two partons in the initial state, giving the 
number of initial-final state clusterings A'lp = 2 n (n + 1). 

- The number of final-final state clusterings is given by n (n + 1). For each of these clusterings, we 
have n — 1 remaining possible recoil partons. For the number of final-state clusterings, we hence 
obtain A'pp = (n — 1) n (n + 1). 

- Accordingly, the total number of possible clusterings is A'tot = -^if + -^ff = n{n + 1)^. 
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- Following these observations, the fraction of final-final state branchings to be generated is then 
given by A^ppZ-^tot = {n — l)/(n + 1). Similarly, the fraction of initial-final state branchings thus 
is A^ip/iVtot = 2/(n + 1). As we see for a very large number of jets, the fraction of pure final-state 
branchings tends towards one as expected. 

Note that in Eq. (I20l ) we see that each dipole term is "averaged" by exactly the number of final-final state 
dipoles, A'pp, or the number of initial-final state dipoles per beam, ^ A'lp. 



B The Explicit Jet Events 

The jet configurations used to calculate the results in Table[T]are given in this appendix, together with some 
kinematic properties of the particular n-jet event. These properties are calculated from the final-state jets. 

• For 71 = 2, we have: 

as Q(mjj)^ = 0.0854525, (mjj) = 2320.1 GeV , 

min(mjj ) ^ 2320.1 GeV , max (mjj) = 2320.1 GeV (35) 

and the momenta (in GeV) are given by 

XaPa = (651.429,651.429,0,0), 
XbPb = (2065.78,-2065.78,0,0), 

Ji = (988.026,4.76957,150.427,976.496), 

J2 = (1729.19,-1419.12,-150.427,-976.496). (36) 

• For ji — 3, we have: 



I = 0.08936, (mjj) = 1546.75 GeV 



1 

min(mjj) = 837.178 GeV, max(mjj) = 2342.36 GeV (37) 



and the momenta (in GeV) are given by 

XaPa = (1002.78,1002.78,0,0), 
xbPb = (1789.36,-1789.36,0,0), 

Ji = (1203.62,-339.322,1151.26,-90.3834), 

J2 = (1243.44,-297.018,-1187.58,-218.141), 

J3 = (345.076,-150.236,36.3206,308.525). (38) 
For 71 ~ 4, we have: 



( o("ijj> = 0.0873363, (mjj) 1899.39 GeV 



1 

min(TOjj) = 906.006 GeV, max(mjj) = 2636.08 GeV (39) 
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and the momenta (in GeV) are given by 

Xa Pa = (2293.54, 2293.54, 0, 0) , 

XbPb = (2359.46,-2359.46,0,0), 

Ji = (1725.71,877.743,982.549,1114.55), 

Ja = (954.182,356.832,-669.803,-578.359), 

Jg = (1037.51,-521.634,-809.127,-386.844), 

Ji = (935.596,-778.871,496.381,-149.35). (40) 



• For n = 5, we have: 

asQ(mjj)^ = 0.0921983, {rrijj) = 1074.2 GeV , 

niin(mjj) = 365.996 GeV , max(mjj) = 1794.95 GeV (41) 

and the momenta (in GeV) are given by 

x^Pa = (883.985,883.985,0,0), 

XbPb = (3263.37,-3263.37,0,0), 

Ji = (684.733,-446.345,-519.146,11.084), 

J2 = (780.483,-90.2618,684.399,-364.149), 

J3 = (1081.8,-949.154,502.292,130.739), 

Ji = (458.187,-330.582,216.105,-232.271), 

J5 = (1142.15,-563.042,-883.651,454.598). (42) 



• For n = 6, we have: 

asQ(mjj)^ = 0.0898744, (mjj) = 1470.28 GeV, 

niin(mjj) = 372.579 GeV , max(mjj) = 2307.33 GeV (43) 
and the momenta (in GeV) are given by 





= (2711.42,2711.42,0,0), 


XbPb 


= (2989.73,-2989.73,0,0), 


Jl 


= (1305.48, -936.767, -465.008, 781.36) , 


J2 


= (416.233, 43.8532, 163.26, -380.359) , 


h 


= (1255.05, 418.867, 211.911, 1163.96) , 


J4 


= (809.408, -755.773, -267.267, -111.885) 


h 


= (1029.43, 606.551, 515.234, -652.956) , 


J% 


= (885.548, 344.959, -158.13, -800.121) . 
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• For n = 8, we have: 



[^2^"^«^j = 0.0956911, {mjj) = 861.375 GeV, 
min(mjj) = 239.955 GeV , max(mjj) = 1483.52 GeV 



(45) 



and the momenta (in GeV) are given by 



•^a Pa 


= (1566.71,1566.71,0,0), 


XbPb 


= (3315.09,-3315.09,0,0), 


Jl 


= (771.638, 72.2583, 756.755, 132.381) , 


J2 


= (799.106, -447.319, -526.283, -401.875) , 


h 


= (362.664, -225.218, 127.162, -254.228) , 




= (455.86, 294.142, -344.111, -53.627) , 




= (336.571, -176.658, 48.6816, 282.315) , 




= (897.978, -705.684, 548.115, 89.1335) , 


Jl 


= (897.601, -802.832, -399.01, 44.0329) , 


Js 


= (360.374, 242.933, -211.31, 161.867) . 



• For n = 10, we have: 

as Q(mjj)^ = 0.0950109, (mj^) = 913.794 GeV , 

niin(mjj) = 237.586 GeV , max{mjj) = 1830.56 GeV (47) 
and the momenta (in GeV) are given by 



Pa 




(2827.46,2827.46,0,0) , 


XbPb 




(3322.41,-3322.41,0,0), 


Jl 




(547.589, 302.187, 339.049, 305.912) , 


J2 




(956.324, -806.585, 389.162, 335.454) , 


Js 




(350.22, -95.7603, 228.577, 247.461) , 


Ja 




(259.829, -26.1444, -35.2505, -256.096) , 


J5 




(314.102, -35.8975, 4.96171, 312.004) , 


J% 




(986.504, -606.912, 589.914, -506.804) , 


Jl 




(889.429, -326.063, -798.032, -218.89) , 


Js 




(891.285, 649.951, -560.226, -241.04) , 


J9 




(518.799, 168.343, -476.506, 117.284) , 


Jio 




(435.784, 281.925, 318.35, -95.2848) . 



• For n = 15, we have: 

as(^{mjj)] = 0.0993949, (wj-j) = 633.545 GcV , 





147.1 GeV, max(mjj) = 1497.81 GeV 



(49) 
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and the momenta (in GeV) are given by 

XaPa = (3403.1,3403.1,0,0), 

XbPb = (3096.06,-3096.06,0,0), 

Ji = (349.237,141.394,-212.337,-238.51), 

Js = (272.807,-73.295,214.113,-152.339), 

Js = (465.556,289.402,279.374,234.392), 

J4 = (269.21,54.8123,-263.571,-0.174956), 

J5 = (451.568,254.001,87.8387,-362.88), 

Je = (437.866,292.297,303.293,-119.594), 

J7 = (368.167,-241.648,223.073,165.504), 

Jg = (355.366,-41.1757,-217.718,-277.828), 

Jg = (328.823,-87.4318,-296.434,-112.281), 

Jio = (458.585,-188.764,-88.852,408.379), 

Jii = (332.789,159.417,-133.945,259.603), 

J12 = (572.199,-470.06,-321.957,-52.8996), 

Ji3 = (746.136,-332.113,245.182,-621.534), 

Ji4 = (830.982,518.077,212.244,614.069), 

Ji5 = (259.874,32.133,-30.3028,256.093). (50) 
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